Influence of the nucleon-nucleon collision geometry on the determination of the 
nuclear modification factor for nucleon-nucleus and nucleus-nucleus collisions 
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The influence of the underlying nucleon-nucleon collision geometry on evaluations of the nuclear 
overlap function (Tab) and number of binary collisions (N co u) is studied. A narrowing of the spatial 
distribution of the hard-partons with large light-cone fraction x in nucleons leads to a downward cor- 
rection for Ncou and Tab , which in turn, results in an upward correction for the nuclear modification 
factor Rab- The size of this correction is estimated for several experimentally motivated nucleon- 
nucleon overlap functions for hard-partons. It is found to be significant in peripheral nucleus- nucleus 
and nucleon-nucleus collisions, and are much larger at the LHC energy of -^/s = 5.5 TeV than for the 
RHIC energy of -^5=200 GeV. The implications for experimental measurements are also discussed. 

PACS numbers: 25.75.-q 



I. INTRODUCTION 

In experiments at the Relativistic Heavy Ion Collider 
(RHIC), modification of hard-scattering processes is usu- 
ally quantified via the nuclear modification factor. For 
collision between nucleus A and B (A-B collision), it is 
defined as 
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where Tab is the nuclear overlap function calculated 
as the convolution of the thickness functions Ta,b(&) = 
J pa,b(&, z)dz for A and B, 



Tab(&) = J dsT A (s)T B {s-b), 



(2) 



and (Tab) is the average nuclear overlap function for the 
corresponding centrality bin 
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In Eqs. [T] and [3J (N co il) is the average number of in- 
elastic nucleon-nucleon (n-n) collisions with inelastic n- 
n cross section At the nominal RHIC energy of 

y/s = 200 GeV, cr™ 6 ' = 42mb = 4.2/m 2 . The estimation 
of N co u, Tab and other geometrical quantities are usu- 
ally obtained either numerically via an optical Glauber 
approach or statistically via a Monte-Carlo Glauber ap- 
proach (see Appendix [Al . More details on the Glauber 
model can be found in 

The definition of Tab m Eq. [U ignores the fact that 
nucleons are extended object and the multiplicity of a 
n-n collision also depends on the impact parameter 6 nn . 
This can be accounted for by introducing a n-n over- 
lap function t(b nn ) with normalization of J db nn t(b nn ) = 



/ db Ti 



t(b r , 



1, which generalizes Eq. [5] to 



Wong's formula 0: 

Tab(&ab) = J db A db B T A (bA)T B (bB)t{bAB - b a + b B ) 
= J dsdb nn Ta(s)Tb(s — Hab + b nn )t(b„„). 

(4) 

In this notation, Eq. [2] is obtained for the special 
case when the n-n overlap function is a delta function: 
tibnn) = S(b nn ). Eq. [4] shows that the nuclear overlap 
function Tab (and thus iV C oii) depends on not only the 
nuclear impact parameter 6ab, but also the n-n impact 
parameter b nn through t(b nn )- In other words, the 6ab 
and b nn are correlated. Because of this correlation, Eq. 3] 
in general is not the same as Eq.[5]at fixed &ab- However 
if the nuclear thickness function varies linearly within the 
length scale of nucleon size, Eq.[5]is a good approximation 
since the integral of the first term of its Taylor expansion 
vanishes due to spherical symmetry of the nucleon, 
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This approximation is rather precise for central A-B col- 
lisions, but breaks down in peripheral collisions. Further- 
more one can show that Eq. [4] (hence Eq. [2]) obeys the 
following sum rule 
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where cr A °g is the total A-B geometrical cross-section. 
This equation implies that the integral of Tab, i-e. for 
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FIG. 1: Schematic illustration of the collision geometry for nucleon-nucleon or n-n (a), nucleon-nucleus or n-A (b) and nucleus- 
nucleus or A-B (c) collisions, a) Collision geometry for peripheral (left) and central (right) n-n collisions. The dark (light) 
region indicates the transverse distribution of hard-partons and soft-partons, respectively (figure adapted from J3|), hence the 
hard-scattering rate is suppressed in peripheral collisions, b) Collision geometry for peripheral n-A collisions. The nucleus 
only overlap with top half of the projectile nucleon, hence the corresponding n-n collisions is bias towards large n-n impact 
parameter, c) Collision geometry for peripheral A-B collisions. The corresponding n-n impact parameter distribution receive 
surface biases from both nucleuses. 



minimum bias events, is independent of the n-n overlap 
function. 

The transverse spatial distribution of the partons in 
a highly lorentz-boosted nucleon depends on the light- 
cone fraction x [!, H| . The distribution of hard-partons 
(x > 10~ 2 ) in the transverse plane (xy-plane) is consid- 
erably narrower than that for soft-partons which define 
the geometrical size of the nucleon (see Fig. la) @. The 
reason is that the transverse profile of soft-partons grows 
with ln(l/x) due to Gribov diffusion @, corroborated by 
the possible tamping of the growth of soft-partons at nu- 
cleon center by saturation effects [7]. As a consequence, 
the n-n impact parameter distribution, i.e. the overlap 
function t(b nn ), is narrower for hard-scattering processes 
than that for minimum-bias n-n collisions. If we denote 
the n-n overlap function for hard-scattering processes as 
ths(b nn ), the corresponding nuclear overlap function and 
equivalent number of n-n collisions for Eq. [1] has a form 
similar to Eq. [J] (see Appendix [BJ : 

t Ab(Pab) = J dsdb nn T A (s)T B (s - b AB + b nn )t hs (b nn ) 

\rhs rjihs inel fc\ 
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The primary goal of this paper is to investigate the sen- 
sitivity of Tab and N co u to different assumptions about 
n-n overlap function for hard-scattering processes. The 
influence of this sensitivity on the interpretation of the 
nuclear modification factor Rab is also discussed. 



II. RESULTS AND DISCUSSION 

To help visualizing the level of correlation between bAB 
and b nn , we plot the b nn distributions for various ranges 
of bAB m Au-Au collisions in Figure O These distri- 
butions are calculated using the MC Glauber approach 



(see Appendix [A"|l . A n-n total inelastic cross section of 
a rm = 42 mb is used in the calculation, which corre- 
sponds to a geometrical radius of r n = \J a^ 1 /tt/2 = 
0.5781/m for nucleons. A n-n collision is considered to 
occur when their distance in the xy-plane is less than 
dmax = 2r„ = 1.156/m. The individual distributions 
in Figure [2] have been re-scaled to match each other at 
6, m =2.5fm. For minimum bias Au-Au events, the distri- 
bution increase linearly with b nn . This is so because all 
possible xy positions of a given nucleon in one nucleus rel- 
ative to the second nucleus are uniformly sampled (con- 
sistent with the sum rule Eq. [5]) . Thus the corresponding 
b nn distribution has the same shape as that for pure n-n 
collisions. By contrast, the distributions for centrality se- 
lected Au-Au events are not linear functions of b nn . For 
peripheral Au-Au collisions, the n-n distribution has a 
concave-like shape, implying that the selection of large 
bAB results in a bias towards peripheral n-n collisions 
(large b nn ). For central Au-Au collisions, the n-n distri- 
bution has a slight convex-like shape, implying that the 
requirement of a small bAB value leads to a bias towards 
slightly central n-n collisions. 

The origin of the bias can be further illustrated via 
peripheral n-A collisions. According to Eq. [4] the n-n 
impact parameter distribution for fixed b n A should be 

f(b nn ,b nA ) = J dsd(j)T B {s)TA{s-bnA + bnn)b nn 

= J d(f> T A (b n A - Kn) Km (7) 

after exchanging subscript A and B and setting T B = 
8{s). When the distance of nucleon from the nuclear 
surface is of the order of nucleon size, i.e. \b n A — R\ < r n 
(R is the nuclear radius and r n = d max l2 is nucleon 
radius), the projectile nucleon see many more nucleons at 
the inner side of the target nucleus than the side close to 
the surface (see Figure [TJd) , leading to the bias shown in 
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FIG. 2: a) Nucleon impact parameter distributions of n-n 
collisions for centrality selected Au-Au events. From top to 
bottom, the curves show several Au-Au impact parameter 
selections, ranging from 0-1 fm to 14-15 fm. The minimum 
bias Au-Au distribution and pure n-n distribution fall on top 
of each other. The two vertical dashed lines at 1.156fm and 
1.514fm indicate the radius of the n-n inelastic cross-section 
at Vs=200 GeV (42mb) and 5.5 TeV (72mb), respectively, 
b) The ratios between various centrality selected Au-Au b nn 
distributions to n-n. 

Figure [2j This bias is larger for A-B collisions, since the 
overlap function Tab is the convolution of two thickness 
functions, which varies more rapidly towards the edge of 
the overlap region. Furthermore, the magnitude of such 
bias should also be sensitive to the nuclear profile. Due to 
the surface diffuseness of the Woods-saxon nuclear profile 
(used in this analysis), the range of b\B affected by this 
bias is expected to be larger than that for hard-sphere 
nuclear profile. 

Measurements of hard inclusive scattering processes 
have provided a detailed picture of the longitudinal mo- 
mentum distribution of partons in a nucleon. Unfortu- 
nately, information on the transverse distribution of par- 
tons in a nucleon, or the impact parameter dependence 
of the hard-scattering process in n-n collisions, is rather 
limited. In general, it is a function of x and Q 2 and can be 
described by projecting the generalized parton distribu- 
tion function (GPD) JEHU- To estimate the magnitude of 
the bias, we tried several different n-n overlap functions. 
The default overlap function is a step function defined as 

tmb{bnn) oc 0(b nn - 2r„), (8) 



which treats all n-n collisions on equal footing in calcu- 
lating the number of minimum bias n-n collisions. It is 
the default approach commonly used in MC Glauber cal- 
culations of Tab and N co u We compare it to the 
following four parameterizations of hard-scattering over- 
lap functions (see Figure [3]) , all of which have a width 
narrower than the default. The first one assumes that 
the hard-scatterings are restricted to a core region which 
is half the size of the n-n overlap, i.e. 

t\ a {b nn ) oc 9(b nn - r n ). (9) 

This overlap profile has no rigorous physics foundation, 
other than the fact that it provides a simple and intuitive 
estimation of the bias. The second one assumes that the 
overlap function is the folding of two Gaussian profiles 
for hard-partons with the width r n . 

1 _" 2 ™ 

thsibnn) « ^-e ^6(b nn -2r n ). (10) 

This is the n-n overlap function used by PYTHIA [§], 
except that we truncate the overlap function at 2r„, so 
the average width is narrower than the default (about 
0.7r„ in this case). The third one assumes that the hard- 
partons are uniformly distributed in a hard-sphere nu- 
cleon with radius r n , the corresponding overlap function 
is 




(11) 

The last one is the dipole formula for gluons spatial dis- 
tribution taken from [3j. It is derived from fits to J/\f' 
photo-production data at HERA and FNAL, 

tUbnn) oc ^| (^)' K 3 (m g b nn ). (12) 

where K3 is the modified Bessel function, and the mass 
parameter m 2 ~ l.lGeV 2 , which depends weakly on Q 2 
and assumed to be constant in this analysis. 

With these input distributions in hand, it is fairly 
straight forward to evaluate the resulting differences in 
^coii anc l N co u using the MC Glauber approach. The 
results obtained for Au-Au collisions are shown in Fig- 
ure |U they are calculated for Eq. [5] with er,™ e( = 42 mb 
(for v / s=200 GeV). The results are plotted against either 
a) Au-Au impact parameter, b) N par t, or c) centrality in 
10% steps (sliced according to Au-Au geometrical cross- 
section). Here, we shall focus the discussion on the mid- 
dle panel. One can see that the N^ H is always smaller 
than the nominal N co u value for small N part . This is 
easily understood since N co u is calculated for a flat prob- 
ability distribution while all others are calculated for a 
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FIG. 3: The shapes of the impact parameter distribution of 
the overlap function for hard-scattering processes at yfs — 
200 GeV; the solid vertical line indicates the range of impact 
parameters for the minimum bias n-n condition. For tje LHC 
energy of yfs=b. 5 TeV, we simply increase the nucleon radius 
r n from 0.578fm to 0.757fm, which stretches all distributions 
(except the dipole formula) horizontally by 30%. 



distribution which peaks at small b nn . The bias is sizable 
at Np ar t < 20 — 50, corresponding to > 60% centrality 
range. The bias for the Gaussian profile is smallest since 
it's t nn is broader than for the other scenarios. For the 
more realistic case of hard-sphere and dipole overlap, the 
bias grows to about 5-15% for N part ~ 10—15 which 
corresponds roughly 70-80% centrality bin. Note that 
the corrections all cross at N part ~ 150 (or bAuAu ~ 8 
fm) and stay slightly above unity for central collisions. 
This is expected since there is no bias for minimum bias 
Au-Au selection, as mandated by the sum rule Eq. [5] 

Based on the trends shown in Figure^ it is clear that 
the bias should grow with an increase of total n-n in- 
elastic cross section. At the LHC energy of 5.5 TeV, the 
total n-n cross section is estimated to be er^" ei w 72mb = 
7.2/m 2 about 75% higher than that at the nomi- 
nal RHIC energy. This corresponds to d max = 1.514/m 
and r n = 0.757 fm. The results obtained for this con- 
dition are summarized by Figure which shows about 
60% larger bias on N co u when compared to that for the 
RHIC energy. 

As discussed before, the same bias should also exist for 
p-A collisions, but to a lesser extent. In Figure[Hl we show 
the results of the calculation for deuteron-Au (d-Au) col- 
lisions at y^s = 5.5 TeV. Indeed the bias is smaller com- 
pared to Au-Au (Figure [5]) at the same N part . However 
the available range of N part covered by d-Au collision is 
small. So the correction factor for given percentile cen- 
trality bin, sliced according to geometrical cross-section, 
turns out to be larger. For the (N part ) values for 60-100% 
centrality range is only about 4, for which the correction 
could be as large as 15%. 

The results in Figure HUH underlines the importance 
of the nucleon-nucleon collision geometry for proper in- 
terpretation of the collision geometry in p-A and A-A 



collisions. Depending on the assumed hard-parton pro- 
file, the bias on the N co u (hence Tab) could become siz- 
able, especially at LHC energy where the total n-n cross- 
section is significantly larger, and in peripheral p-A and 
A-A collisions where the correlation between bAB and 
b nn is important. In such cases, the full formula for Tab 
(Eq. 3]) and proper n-n overlap function for hard-partons 
are needed for proper evaluation of the nuclear modifica- 
tion factor. The magnitude of the bias is usually smaller 
than the typical systematic error quoted by the RHIC ex- 
periments |l| . Also the experimental triggering efficiency 
for peripheral A-A and p-A collision is low, for exam- 
ple, the PHENIX experiment triggers on 92% of Au-Au 
collisions [HI and 88% of the d-Au collisions [l| , hence 
the relevant correction factor for experimentally accessi- 
ble centrality range could be 30-50% smaller than what 
has been calculated here. We should point out that a 
related bias has been considered before by the PHENIX 
Collaboration [l[, where a Gaussian profile was used as 
part of the estimation of the systematic errors on N co u . 
However that study was not motivated by the narrowing 
of hard-parton spatial distribution, rather it was used to 
account for uncertainty of the nucleon matter distribu- 
tion which are dominated by soft partons. 

In summary, we discussed a bias in the calculation of 
the N co u and Tab , and its influences on the nuclear mod- 
ification factor Rab within the Glauber formulism. The 
bias is caused by the difference in the spatial profile be- 
tween the hard-partons which contribute to the high px 
yield, and soft-partons which determines the n-n inelastic 
cross-section. The much narrower spatial profile of the 
hard-partons biases the N co u and Tab to values smaller 
than those obtained using the uniform n-n overlap func- 
tion. The magnitude of the bias is sensitive to width of 
the overlap function for hard-partons. A crude estima- 
tion for Au-Au collisions at RHIC energy leads to about 
5% (assuming a Gaussian overlap function)-15% (assum- 
ing a hard-sphere overlap function) downward corrections 
for 70-80% centrality bin and becomes significantly larger 
for more peripheral collisions, which may account for part 
of the suppression seen for very peripheral Au-Au tt° 
data For Npart > 150, the bias is positive, but is 

less than a few percent. Since the bias increases with the 
total n-n inelastic cross-section used in the calculation. 
We estimate the bias could be 60% larger at LHC energy 
compare to RHIC energy, and has to be properly taken 
into account in Glauber calculation for both Pb-Pb and 
p-Pb collisions. 

The author wishes to thank Roy Lacey, Michael Tan- 
nenbaum and Cheuk-Yin Wong for valuable discussions. 
This research is supported by NSF under award number 
PHY-0701487. 



APPENDIX A: OPTICAL VS MC APPROACH 

The calculations of Glauber variables, including (Tab) 
and (iV co ii), are carried out either numerically using a set 
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FIG. 4: The ratio of N CI ^ U to N co ii calculated for Au-Au collisions for the four n-n overlap functions described by Eq. [9lTl2l 
plotted as function of a) impact parameter bAuAu, b) N par t and c) centrality bins in 10% step, assuming n-n inelastic cross 
section of 42mb (for RHIC energy). 




FIG. 5: The ratio of N C< ^ H to N co u calculated for Au-Au collisions for the four n-n overlap functions described by Eq. [9lTl2l 
plotted as function of a) impact parameter bAuAu, b) N par t and c) centrality bins in 10% steps, assuming n-n inelastic cross 
section of 72mb (for LHC energy). 




FIG. 6: The ratio of N^ u to N co ii calculated for deuteron-gold (d-Au) collisions for the four n-n overlap functions described by 
Eq . I51TT21 plotted as function of a) impact parameter bAuAu , b) N par t and c) centrality bins in 10% steps, assuming n-n inelastic 
cross section of 72mb (for LHC energy). The Hulthen wave function for deuteron is used. 
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of equations like Eq. [2j[4] ( "optical approach" ) or statis- 
tically via Monte-Carlo method ("MC approach"). The 
MC approach is what has been used by RHIC experi- 
mentalist. Comparing with Optical approach, it takes 
into account position fluctuations of nuclcons in the nu- 
cleus, which is shown to be important in peripheral colli- 
sions. In the MC approach, A-B collisions are generated 
randomly in the xy-plane. Within each event, nucleus 
A and B are populated randomly with nucleons accord- 
ing to the Woods-Saxon nuclear profile. All possible n-n 
combination between the two nucleus are considered. In 
the simplest version at RHIC energy, a n-n collision is 
considered to happen when the ir distance in the xy-plane 
is less than d ma x — \J <J'mf / ir = 1.156/m (hard-sphere 
assumption, corresponding to t(b nn ) = 6{\b nn \ — d rnax )). 
The number of n-n collisions in a given event is 



A^coll 



E 



t(\n-?j\) 



(Al) 



The Tab is then calculated as T AB = N coU /a 1 ^ 1 . The 
(Aeon) and (Tab) are obtained by averaging over many 
events falling in a given centrality definition. More details 
can be found in 



a nn T AB( b AB) 



(B2) 



where f(b nn ,bAB) is n-n collision density in impact pa- 
rameter space defined by Eq. [71 and 



t ab = J dsdbnn T A {s)T B (s - b AB + b rm )t hs (b nn )(B3) 

has the same form as Eq.4 except that the overlap func- 
tion is replaced with the one for hard-scattering process. 

The n-n hard-scattering rate is simply the hard- 
scattering profile averaged over all n-n events in vacuum, 



ryhsx = 
\ nn I 



Vnnf db nn 2wb nn t hs (b nn ) _ a, 
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where the normalization J db nn 2irb nn th s (b nn ) = 1 is 
used. 

The correct scaling factor should be the ratio of the 
hard-scattering rate in centrality selected A-B collisions 
to that of the minimum bias n-n collisions: 
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APPENDIX B: SCALE FACTOR FOR 
HARD-SCATTERING PROCESS 

In this section, we derive the expression for correct 
scaling factor, N^ a , for hard-scattering process. The av- 
erage hard-scattering yield in A-B event at a fixed impact 
parameter can be expressed as 



Yab = ( E M&-6il)<# > £ 

i&A,jeB 



(Bl) 

a nn j dbnn f(bnn,bAB)ths(bnn) 

(T^n / dsdb nn T A (s)T B (s - b AB + b nn )t hs (b nn ). 



This equation is similar to the original definition Eq. [31 
except that a different n-n overlap function is used. 

A few words about N£* u are in order. First, according 
to Eq. [51 there is no bias for minimum bias A-B event 

selection, i.e. (Nj£ H ) = (N coU ) = AB J B W- . Secondly, 
the n-n hard-parton overlap function th s depends on the 
x and Q 2 , thus in general, N^ u could depend on the 
Pt- Last but not the least, the magnitude of the hard- 
scattering cross-section <j^ n in Eq. IBll is not important, 
what is the important is it's spatial distribution. In fact, 
(Tnn cou ld be arbitrarily small, and yet would still cancel 
between A-B and n-n (Eq. IB5[) . For our MC Glauber 
calculation, a convenient choice is to choose a^ n — a^ 1 , 
then Eq. IB1I would directly give N^ c 
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